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^ '. Abstract 

A model is developed for electromagnetic form factor of the pion. One-loop cor- 

1 

^ ■ rections are included in the linear u— model. The /)— meson contribution is added in 
an extended VMD model. The form factor, calculated without fitting parameters, 

I is in a good agreement with experiment for space-like and time-like photon mo- 

I menta. Loop corrections to the two-pion hadronic contribution aj!^"''^''^^ to the muon 

Cp ■ anomalous magnetic moment are calculated. The optimal value of the u— meson 

. mass appears to be close to the /?— meson mass. 



Oh! 

?H ; 1 Introduction 



X 



It has recently been understood that the pion electromagnetic (EM) form factor is a very 
important physical quantity that plays a key role in the test of the Standard Model at 



ci I the electroweak precision level. The reason is that the production cross section 

na^ ( Ami ^ 



^(e+e-^7r+7r-) = ^(^l-^J , (1) 

where s is the squared total energy in center of mass system, a is the fine-structure 
constant and is the pion mass, at low energies dominates over the other hadronic 
channels and accounts for more than 70 per cent of the total hadronic contribution to the 
muon anomalous magnetic moment (AMM) = {g^ — 2)/2. The recent measurement of 
from Brookhaven E821 experiment [1] has boosted the interest in a renewed theoretical 
calculation of this quantity [2]. 

The main ingredient of the theoretical prediction of a^, which is responsible for the 
bulk of the theoretical error, is the hadronic contribution to the vacuum polarization. 
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The contribution of the tt+tt channel to the electron-positron annihilation process can 
be written in terms of the form factor F„{s) via the dispersion integral [3] 



^iHaa,.) ^ ^(^) 1-^ \F^{s)\'<ls, (2) 

4 Jiml \ S J 



oTT^s Jo 2; + (1 — xjs/m'jj^ 

where m^^^ is the muon mass. 

Conventional strategy of the model-independent evaluation of this integral consists 
in utilization of precise experimental data (at least at low energies, where perturbative 
QCD cannot be reliably applied). However, the announced accTiracy, which is to be 
reached soon in E821 experiment, requires calculation of EM radiative corrections to 
the cross section (1) [4]. Apart from the pure tt+tt" events, EM radiative corrections 
include the 7r"'"7r~7 process, where photon is radiated from the final pions. In the current 
experiments at $ and B factories, based on radiative return method [5], this contribution 
cannot be extracted in a model-independent way'^ and corresponding procedure requires 
model-dependent approaches. This, in turn, stimulates development and study of difi'erent 
theoretical models of pion-photon interaction. The simplest one is the point-like scalar 
quantum electrodynamics (sQED) [7] joined with standard Vector Meson Dominance 
(VMD) model (see, for example, [8]) for description of the 7* vr+vr" transition form 
factor in the p-resonance region. Such a model was used in Ref.[4] for construction of the 
Monte Carlo event generator. 

In the present paper we consider modification of the pion EM form factor in the linear 
(7— model [9] with spontaneously broken chiral symmetry, which includes the nucleon 
sector. The p— meson contribution is added following Refs. [10, 11]. In particular, the p 
coupling to the pion and nucleon is introduced through gauge-covariant derivatives, while 
direct 7p coupling has explicitly gauge-invariant form. We calculate the pion form factor 
in one-loop approximation in the strong interaction and compare with the precise 

data obtained from elastic e~Tr~^ scattering and e'^e~ annihilation in the pion pair. 

We pay attention to the effect of the loop corrections in a^i!i^"''^'^\ In general, as La- 
grangian of the (j— model contains Lagrangian of sQED as a constituent part, one can 
expect that the difference between predictions of a— model-|-VMD and sQED+YMD 
is small. Indeed, it follows from our calculation that the loop corrections increase the 
low-energy part of the right-hand-side of eq.(2) by about 2 per cent, as compared with 
sQED+YMD. 



2 Formalism 
2 . 1 Lagrangian 

Lagrangian of the model consists of the two parts: C — C^^^ -\- C^'^^ . The first one is 
Lagrangian of the chiral linear cr— model [9] with an explicit symmetry-breaking term c0. 

^Even in direct scanning experiments a model-independent treatment of the 7r+7r~7 events suggested 
in [6] seems too complicated to be used in the near future 
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After spontaneous breaking of chiral symmetry and re-definition of the scalar field via 
4> = a + V, where v = (0) is vacuum expectation value, Lagrangian takes the form 



= N(i^-rn^)N+^-[(daf-rny]+l[(dnf-rnln'] 



+ const, (3) 



where N, tt and a are the fields of the nucleon, pion and meson with vacuum quantum 
numbers, respectively, is the coupling constant, A is parameter of the meson potential, 
and ^ = d'^Jn, (da)'^ — d'^adfj^a, etc. All parameters of the model are related via 

mN^gnV, ml^2Xv^ + ml, ml^^. (4) 

Moreover, in the tree-level approximation v = where /jr = 93.2 MeV is the pion weak- 
decay constant. More details on the cr— model can be found, for example, in [12], Ch.5, 
sect. 2. 6. 

The second part of Lagrangian includes interaction with the EM field A^, and the 
field of the neutral p— meson. This interaction can be obtained by using the minimal 
substitutions 

^Mtt" ^ ^'^tt" + (eA'' + (a, 6 =1,2, 3), (5) 

where e is the proton charge, gp is the coupling constant, and is the third component of 
the Pauli vector r = (ti, T2, T3). In addition we include the direct coupling of the photon 
to the p— meson in the version of VMD model from Refs. [10, 11]. In this way we obtain 

- {eA^ + gpp^) (tt x di^i:)^ + (e^^ + gpP^.fi^'' - ttJ) 
-g^NY^-^Npp - eNY^^NA, pp,F^\ (6) 

Here = d^A^ - d^A^, p^^ = d^p^ - d^p^ and fp determines the 7p coupling. In eq.(6) 
we assumed equal coupling constants of the p to the pion and the nucleon in accordance 
with the universality hypothesis of Sakurai (see, e.g., [12], Ch.5, sect. 4). At the same 
time the 7p coupling constant fp docs not necessarily coincide with gp. Lagrangian (6) 
obeys gauge invariance because of the form of the 7p coupling. We should mention that 
the nucleon contribution is also included in Lagrangian (6), contrary to [10, 11]. 



2.2 Count erterms and renormalization 

Since one of the purposes of the present paper is to take into account loop corrections to 
the pion EM vertex one needs to specify the way of renormalization of the parameters. We 

use the conventional approach and assume that in Lagrangians (3) and (6) enter the "bare" 
fields, coupling constants and masses, to be marked by the subscript "0" . The bare fields 
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require rescaling: (TTo.ao) = ^fZ^ijx a) Nq = ^JZls|N, = y^Zpp^ and Aq = ^/Z^A^, 
where Z^^^Z^, Zp and Za are the respective wave-function renormahzation constants for 
the pion (or sigma), nucleon, rho and photon. The procedure for obtaining Lagrangian 
of counterterms is known (see, for example, [13], Ch.lO). For C^^^ the corresponding 
counterterm Lagrangian reads 

= Sz^m^N - S,^vNN - 6,^N {a + t^,T7r) N + ^6z^ [{dTrf + {daf] 

-\{5p,^bxv^)v-5^a + co-n&X. (7) 

There are six constants 5^^, 5^^, 5^, 5^^, which can be fixed by imposing, in general, 
six conditions on the Green functions. In the calculation of the pion EM vertex only one 
constant hz^ will be needed (see sect. 2.3). 

For Lagrangian (6) one can define first the physical values of the electric charge e = 
eQ^fZA and the rho coupling gp = g^jQyfZ'p. It is also convenient to introduce ihp — 
\f^pf^pQ (p— meson mass in the absence of the coupling to pions). The counterterm 
Lagrangian can be written as 

4? = -^z}-F,.F^^-6z\p,.P^'' 

-5z^ {eA^ + QpP^) (7f X a^7f)3 + 5zA^.A^ + ^pPm)'(^' - ^3) 
-5z.9pNY^-^Np, - 5z,eN^^^^NA, -5^^- p,.F^r (8) 

It is seen, that in general one needs three additional constants Sz^i ^Zp and 5/^, once 5z^ 
and 5zf^ are fixed. Finally, the total Lagrangian is the sum^ 

£ = (9) 



2.3 Contribution to pion EM form factor from the sector of 
(7— model 

Feynman rules for Lagrangian (9) are obtained according to the standard prescriptions 
[14]. The counterterm constants can be found by imposing the following constraints on 
the self-energy operators of the pion, sigma-meson and nucleon respectively 



S7r(rr4) 



d 

dp 



^aiml) = 0, 



(10) 



lp=mN 



= 0. 



These conditions imply that the pole positions of the pion, nucleon and sigma propagators 
are located at the physical mass of the pion, nucleon and sigma respectively. In addition, 
the residue of the pion and nucleon propagators is unity, ensuring the absence of the 
renormahzation for the external pion and nucleon (but not for the external sigma-meson). 

*The mass in is replaced by nip. 
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We also impose the condition (a) = 0, which is ensured by requiring that the so-called 
tadpole diagrams vanish. Correspondingly, the tadpole diagrams will not contribute to 
quantities calculated below. 

In calculation of the loop integrals we used the dimensional regularization method 
(see, e.g., [13], Appendix A.4). Exploiting the conditions (10) we find the constant Sz^- 



47r2 



/ _ In + mlx{l - x) _^ (ml - mlfx{l - x) 



A2 



dx, 



Atvtt = - "7,^2^(1 - x), 

2 



A^ra = mix + ml{l - xf. 



(11) 



(12) 



In these equations /g = 7^ + In 47r, e = 4 — D — > 0, where D is the space-time 

e 

dimension, 7^ 0.5772 is the Euler constant and A is an arbitrary scale mass, which 
drops out in the physical observables. 

The one-loop contributions to the pion EM vertex coming from the o"— model are shown 
in Fig. 1. Using the isospin structure of the vertices, or negative charge-conjugation parity 
of the photon, one can show that the diagrams "e,f" and "g" vanish. The counterterm 
"h" cancels divergences coming from the loop contributions "b" and "c" , while the con- 
tribution "d" is finite. 

In general case of the off-mass-shell pions, the EM vertex r^^(pi,p2, ?) for the process 
7*(g) 7r"(pi) -|- 7r''(p2) has the form 

-iKbiPuP2.q)^e'^ [F{plpW){p,-p,Y-rG{plplq'){p,-rPir] , (13) 

with scalar functions F{pl,pl, q^) and G{p\,p\, q^). On the mass shell, Pi = P2 = ^j^, the 
function G{m^,m^,q^) drops out, while F((m'^,rn^,q^) becomes equal to the pion form 
factor F^(g2). Denoting the loop corrections by i^F'^"'\q^) we find 



FW(g^) = l + AFW(a + 5^^, 
where the total correction is finite and equal to 

m\ — mlx{l — x) 



(14) 



AF^'^\q')+Sz^ 



47r2 
+ 



In 



mif-q^x{l-x) ^A^n 
^ y(l-y)(m2-0\ 



Jo ^ 



{ml - mlf 
AmlrA 



''N'-^ncr 



]) 



dx 



y^ml 



ANn{y,x) Am%A^^{y,x) J 



AnAv^ x) = - q y x{1 - x) - m^y{l - y), 



A^a{y, x) = ml{l -y)- y^{q^x{l - x) 
and Ajvtt, A^^ are defined in eq.(12). 



(15) 
(16) 



2.4 Contribution to pion EM form factor from the p— meson 



The contribution to the pion EM form factor from the p— meson can be written in the 
compact form 



F(^\q') 



9piq 



q^ 



f,{q') q'-rfij- U,{q^) 



+ AFl^P^\q^). 



(17) 



Figure 1: One-loop diagrams contributing to the pion EM form factor in the cr— model. 
Dashed lines depict pion, dotted lines - sigma, solid lines - nucleon, and wavy lines - 
photon. Small crossed circle denotes the counterterm. Diagram "a" corresponds to pion 
form factor in sQED. 



This expression includes several effects coming from the loop corrections which are shown 
in Fig. 2. 

i) The g^— dependent vertex gp{q^) describes loop corrections to the pim coupling, 
which originate from the a— model (Fig. 2, diagrams "a"). These corrections have not 
been included in Ref.[ll]. We can write 

rf)=^?p[l + AF('^)(g2) + 5zJ. (18) 

It is seen that the expression in the square brackets is the same as in eq.(14) and is finite. 
From eq.(18) one obtains 

( 2^ ( 2^ l + AFW(g^) + ^z. 

+ + (19) 

in terms of the constant gp{'m^p). From the experimental width of the p — »• tttt decay, 
Fp^^TT = 150.7 MeV [15], we get \gp{Tn?p)\ = 6.05. To find the real and imaginary part of 
gp{rn?p) one can make use of the relations 

Re gp{ml) = -^=i==|^p(mj)|, Im ^p(m2) = -^=^==|c,^(mj)|, (20) 

A = ImAF(")(m2)[l + ReAF^''\ml) + 6zX^. 

ii) The p— meson self-energy has the structure n^''(g) = {g^^ — q^q^ / q^)Tlp{q^) and cor- 
responds to the diagrams shown in Fig. 2 "b" . It leads to the following exact propagator 
of the rho-meson (Fig. 2, diagrams "c") 

^7(1) = 2 -2 n?2^ + (21) 
^ q'^ — mj — Ilp(g^) g^m^ 

Calculation of the loop integrals in Fig. 2 "b" results in 

np(g^) = q'Hq')-6z,], uj{q') = ^^[-h + I{q% (22) 

/(g') = 2^ [2xln^ + (l-2a;)ln^](l-x)da;, (23) 



b , 




Figure 2: Diagrams "a" show vr, a, loop corrections to pixn vertex; diagrams "b" - 7r,N 
loop corrections to self-energy of p— meson; diagrams "c" graphically represent equation 
for exact p— meson propagator; diagrams "d" -tt,N loop corrections to 7p vertex. Double- 
wavy lines depict p— meson. The corresponding counterterms are indicated by crossed 
circles. 



with Atv = fnj^ — — x) and = — g^x(l — x). The self-energy has the 

logarithmic divergence and needs renormalization. The authors of Ref. [11] renormalized 
the self-energy by applying a dispersion relation with two subtractions. We prefer an 
alternative method of counterterms which is expressed in eq.(22) through the constant 
6zp- One can fix the latter from the constraint on the self-energy at the physical mass rrip 
of the p— meson: 

■^ReUp{q%2=mi = (24) 
5z^ = Zp-l = Re iu{ml) + mjRe ^'(mj), (25) 

where Re u}'{q'^) = ^Re ^{q'^)- It is seen from eqs.(22) and (25) that the self-energy 
np(g^) = q'^[u!{q'^) — Re t^(m^) — m^Re u'^m"^^)] is finite. Near the physical mass the latter 
has the expansion 

Re n,(g2) = -rr^jRe a;'(mj) + 0{{q' - mj)^), (26) 

and therefore the coupling Qp is not renormalized due to the self-energy loops [11]. There 
is also a finite mass shift 

— = — m^Re ^^'(m^). (27) 

For definition of rhp see the paragraph before eq.(8). 

Above the two-pion threshold the self-energy acquires an imaginary part coming from 
the pion loop (the third diagram in Fig. 2 "b" ) . Namely, at q^ < 4m^ the imaginary part 
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and the g^-dependent p — > tttt decay width read respectively 

487r ^ g2 

r,{q') = -Imn{q')/y^. (28) 

iii) Closely related to the self-energy are the loop corrections to the jp coupling con- 
stant, shown in Fig. 2 "d". The sum of all contributions is proportional to the tensor 
g^"^ — q^q'^/q'^ similarly to the tree-level term. Introducing the g^-dependent vertex one 
obtains 

where can be fixed by requiring that on the mass shell, q^ = m^, the coupling /p(m^) 
is related to the p e^e~ decay width. The experimental width Tp^^+e- — 6.77 KeV 
[15] is reproduced with |/p(m^)| fa 5.03. Prom eq.(29) we find 

+ -[u;{ml)-u{q% (30) 



fp{.<f) fpiml) Qp 

Here the real part of the constant fpinr^j) is determined from \ fp{m'j^)\ and the imaginary 
part, Im fpirrij^) = |/p(m^)plm uj{m'^p) / Qp. It is seen from (30) that the effective 7p vertex 
is finite. A similar procedure for this vertex was used in [11], although only the real part 
of fpim^p) was taken from experiment. 

Let us also mention that in calculation of Ilp{(f) and fp{(f) we used |g'p(m^)| (instead 
of Qp) in order to get the correct width of the p— meson. 

iv) The last term in eq.(17) describes the p — u interference due to the EM effects [8]. 
The explicit form of the contribution to the pion form factor can be taken from Ref.[16]: 

^Fi^-\q') = -Spjf- . ^ , (31) 



ml, — m?p 



i[mj:^ - mpVp{q^)] ' 



where F^^ = 8.43 MeV is the full decay width of the c^;— meson with mass m^^, — 17.05 
is the 7CJ coupling constant which is fixed from the uj — * e^e^ decay width r^^_»e+e- = 0.6 
KeV [15], and flp^^ ^ —3.8 ■ 10^'^ GeV^ is the mixed p — u) self-energy. 



3 Results and discussion 

Let us first specify parameters of the model. The constant g.,r is determined from the 
trec-lcvcl Goldberger-Treiman relation g.^ = mAr//^ (first equation in (4)), while gp{m^) 
and fp{iTi'^p) are fixed from experiment, as described in sec. 2.4. The a mass is chosen 
equal to the mass of the p, i.e. = rUp, in line with Ref.[17], where a and p are assumed 
to be degenerate. Furthermore, nip = 768.5 MeV and m^; = 782.57 MeV [15]. 

It is interesting to note that calculation of the self-energy of the p— meson gives rhp — 
795 MeV. This value is rather close to the physical mass rrip. In this connection note that 
the authors of [11] fitted Trip from the tttt scattering and obtained 810 MeV. The difference 
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in the above values of rhp is partially due to our taking into account the nucleon loop, 
which was not considered in [11]. 

As mentioned in sec. 2.4, both the real and imaginary parts of the coupling con- 
stants g'p(mp) and fpirri^) have been included. The calculation yields: Re gp{m'^p) — 
6.036, Im gp{ml) = 0.405, and Rc /^(m^) = 4.96, Im /^(m^) = -0.82. Taking into 
account the imaginary parts leads to a small correction to results obtained when setting 
Im gp{ml) = Im /^(m^) = 0. 

Our main results are demonstrated in Fig. 3 and Tables 1 and 2. The calculated 
pion form factor |F7r(g^)p for space-hke and time-hke values of is presented in Fig. 3. 
Apparently the agreement with the data [18] from elastic electron-pion scattering, and 
the data [19] from e^e^ annihilation in two pions is quite good. We emphasize that in 
our approach there are no fitting or tuning parameters. 



There is a strong interference of the two contributions, F^^^ and F^^' , in the total form 
factor. We show in Fig. 3 separately contribution F^"^^ (marked "cr— model"). It is seen 
from eqs.(17) and (18) that the p contribution also includes vr, a, N loops coming from the 
(7— model. Switching off these corrections, i.e. putting AF^'^^q'^) + Sz^ — 0, we obtain 



p{sQED+VMD)(^^2^ ^ _ 



9p 



ml 



+ AFiP'^\q^). 



(33) 



Here the first term corresponds to sQED, the second and the third terms are the p— meson 
contributions. Note, that eq.(33) corresponds to the extended version of VMD of Ref.[ll]; 

in the "standard" VMD model [8] one has the dependence m^/[m^ — — i'm-p^piQ'^)]- 

Our calculation shows that the difference between the form factor, calculated in 
a-model+VMD (eqs.(14) and (17)), and that in sQED+VMD (eq.(33)) is small, and 
therefore the results for sQED+VMD are not plotted in Fig. 3. Nevertheless, the differ- 
ence may show up in the integrated quantity a^'*"^''^^ for the muon AMM. 
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Figure 3: Pion EM form factor for space-like q^ (left) and time-like q^ (right). Experi- 
mental data are from [18] and [19] respectively. 



The calculated values of a}* " are shown in Table 1. In general, the loop corrections in 
the (T— model are important (compare the 1st and the 2nd columns). Their role is however 
diminished in the full calculation, which includes the dominant p— meson contribution 
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sQED 


cr— model 


sQED 
+ VMD 


0"— model 
+ VMD 


cr— model 
+ VMD (/, = g,) 


Ref.[2] 


525 


753 


4667 


4763 


4745 


4774 ±51 



Table 1 : Two-pion contribution a}t " ''^ to the muon anomalous magnetic moment (in 
units 10~^^). The upper integration hmit in eq.(2) is 0.8 GeV^. 



rricj, GeV 


0.4 
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1.1 


1.2 


^(ftad,.)^ 10-11 


4546 


4583 


4640 


4710 


4788 


4867 


4946 


5024 


5099 



Table 2: Dependence of " on the mass of u— meson. 

(the 3rd and the 4th columns in Table 1). The difference between cr— model+VMD and 
sQED+YMD calculations is about 2%. 

In this connection our result can be used to estimate the size of radiative corrections 
due to final-state-photon radiation in the e+e~ — > 7r+7r~7 process. The corresponding con- 
tribution to the muon AMM, af^'^'''''^'^\ is calculated in [4] in framework of sQED+VMD. 
In general, a{e'^e~ v:^t\~^) j a{e^e~ — >• tt+tt-) is of the order a. We expect, that the 
model dependence of similar to that of a'l^"''^'^\ Therefore, the deviation of 

^{had,-K'y) (^^^^(^^jg^^g^ sQED+WMD, from a^J^"'^''^'^'' in a more realistic model, for example, 
in (7— model-|-VMD, is about 2%. The overall model-dependence effect in the contribution 
^{had,TV'y) muou AMM is of the order a x 2% fa 0.015% and is therefore neghgible. 

In the 5th column we show result obtained if wc put /p = g^y in Lagrangian (6). This 
approximation corresponds to the full universality of Sakurai. In this case the renormal- 
ization procedure for the 7p vertex changes and 5/^ in eq.(8) is equal to Szp/gp- The 

numerical results for the form factor and a^'*"'^''^'*, however, change very little, e.g., for the 
integral by about 0.4%. 

As mentioned above, in the calculation we chose the mass mg. = rUp for the cr— meson. 
This particle is associated with /o (400 — 1200)— meson in [15]. In view of its undetermined 
status we study in Table 2 the dependence of the calculated integral a^^"''^'^'' on rricr. As 
can be seen from Table 2, the integral varies considerably. Let us take the value 4774 ±51 
from [2] as a very accurate fit to the experimental integral. Then, for the indicated error 
bars we obtain the mass of the a in the interval from 720 to 850 MeV. The central value 
785 McV is surprisingly close to the p— meson mass. Therefore our calculation agrees with 
the hypothesis of Ref. [17] about degeneracy of a— and p— mesons. 

In the calculations we took into account diagrams, where the p— meson enters only 
on the tree level. In particular, the loops with intermediate p— meson for the 'y7r~^7r~ 
and p7r"'"7r~ vertices are left out. Such contributions can be consistently considered in 
the models, in which p— meson is included together with its chiral partner, the axial- 
vector fli— meson, for example, in the so-called gauged cr- model [20] or chiral quantum 
hadrodynamics [21, 22]. This work is planned for future. 
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4 Conclusions 



We developed a model for EM vertex of the pion. The model is based on the linear 
(T— model, which generates the loops with intermediate pion, sigma and nucleon. The 
/9— meson is included in line with an extended VMD model [11]. The coupling of the p to 
the pion and nucleon is introduced through gauge-covariant derivatives, and the direct 7p 
coupling has gauge-invariant form. The p— meson self-energy and modified 7p vertex are 
generated by the pion and nucleon loops. The renormalization is consistently performed 
using the method of counterterms without cut-off parameters. 

The pion EM form factor, calculated in one-loop approximation in the strong inter- 
action, is in good agreement with the precise data obtained from elastic e~7r+ scattering 
and e^e~ annihilation into 7r"'"7r~. The effect of the cr— model loops turns out to be small. 

We calculated the contribution of the e+e~ — > tt+tt" process to the muon AMM 
moment, a|f'"'^''^\ The calculation agrees quite well (by 0.15%) with recent very accurate 
fit of Ref. [2]. The contribution of the cr— model loops to a^f'"^''^'' is about 2%. 

We also estimated the size of the model-dependence effects in a^i!^"''^''^"'\ the contribu- 
tion to the muon AMM from final-state-photon radiation in e^e~ — > 7r"'"7r~7 process. It 
is about a x 2% 0.015% and is therefore neghgible. Thus our calculation does not 
contradict to the conclusion of Ref. [4] that final-state radiative process e+e~ — > 7r+7r~7 
can be evaluated in scalar QED supplemented with VMD model. 

The only free parameter of the model, which is not fixed from the experiment, is the 
cr— meson mass. Comparison with the fit of Ref. [2] strongly points to the value of this 
mass close to the mass of the p— meson. This conclusion is consistent with Ref. [17], where 
the mesons a and p are assumed to be degenerate. 
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